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Generally,  the  irregular  behaviour  of  a function  f(x)  is  described  by  its  Holder 
exponent  h(x ) quantifying  the  strength  of  the  singularity  at  a point  x.  In  most 
cases,  where  the  function  does  not  contain  oscillating  singularities,  the  Wavelet 
Transform  Modulus  Maxima  method  (WTMM)  allows  a reliable  estimation1,2  of 
the  singularity  spectrum  D(h),  i.e.  the  Hausdorff  dimension  of  the  set  of  all  points 
x with  the  same  Holder  exponent  h(x)  — h.  Methods  based  on  the  WTMM  should 
be  numerically  more  stable  than  direct  methods  such  as  the  Structure  Function 
method1,  since  they  involve  only  weighed-averaged  quantities  and  not  averages 
of  function  increments  and  since  small  errors  are  relatively  less  important  if  the 
calculation  is  restricted  to  maxima. 

However,  in  the  presence  of  oscillating  singularities  the  standard  WTMM 
method  gives  irrelevant  information  on  the  Holder  regularity  of  the  function.  In 
general,  two  exponents  h,/3  are  necessary  to  describe  the  singular  behaviour  of  a 
function  f(x),  namely  the  Holder  exponent  h and  the  oscillation  exponent  fi  de- 
scribing the  local  power  law  divergence  of  the  instantaneous  frequency3,4.  If  f(x) 
contains  oscillating  singularities  the  regularity  of  the  primitive  of  f(x)  depends  on 
/3.  In  this  case  the  Holder  exponent  does  not  increase  by  1 as  in  the  case  of  a 
cusp  singularity  but  by  p + 1.  Thus,  the  singularity  spectrum  of  general  functions 
depends  on  both  exponents,  D{h,(3). 

In  order  to  extract  Holder  exponents  and  to  quantify  at  the  same  time  the 
oscillating  behaviour  we  propose  to  use  a family  of  complex  progressive  wavelets5 

{i’n} 


ip0{x)  = - e u’°/2)e  *2/2,  ipn(x)  = -^n-i(z)  (neZ)  (1) 

ax 

with  an  increasing  number  of  vanishing  moments 


+oo 

J xkipn(x)dx  = 0 for  0 < k < n (2) 

— OO 


For  a sufficiently  large  frequency  wq  (wy  > 5 for  practical  purposes),  ?/;o  can  be 
replaced  by  the  Morlet  wavelet. 

Using  these  wavelets  we  follow  similar  ideas  as  in  an  earlier  paper2,  where  an 
algorithm  was  developed  for  a family  of  real  valued  wavelets  based  on  derivatives  of 
the  Gaussian  distribution,  which  allows  a direct  tracing  of  the  skeleton  of  wavelet 
transform  maxima  lines2.  Using  the  following  abbreviation  for  the  continuous 
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wavelet  transform 
W 


+oo  

/ : = wnf(a,b)  = - [ f(x)ipn(- — - )dx  (a,  6eR,  o > 0) 

a J a 


(3) 


where  a denotes  the  scale,  b the  shift  parameter,  and  %jin  the  complex  conjugate  of 
ipn,  we  set  up  a partial  differential  equation  for  the  wavelet  transform 


d_ 

da 


ioj0~  + ~)  Wnf  = 0 
ob  a / 


(4) 


Describing  the  maxima  line  in  a parametric  form  {a(t),b(t)},  the  motion  along  this 
line  is  given  by  two  ordinary  differential  equations 


with  a constant  C,  which  can  be  integrated  numerically.  Similar  equations  can  be 
derived  for  the  ridges  5.  The  skeleton  of  maxima  lines  can  be  continuously  traced 
up  to  the  desired  accuracy  with  a reduced  computational  effort.  Based  on  these 
maxima  lines  a partition  function  Z{p,q,a)  is  defined 

Z(p,q,a)  = 'Y  (sup  \Wnf(a,bi)\)gdi(a)p  (6) 

a/<a 

zGmax. lines 

where  di(a)  denotes  the  distance  between  the  i-th  and  the  ( i + l)-th  maxima  line. 
In  comparison  with  classical  methods,  the  main  advantage  of  the  proposed  direct 
tracing  algorithm  is,  that  the  supremum  along  each  maxima  line  can  be  easily 
evaluated. 

Exploiting  the  scaling  behaviour  of  the  partition  function  Z(p,q,a)  ~ aT(p'9^ 
and  applying  the  Legendre  transform  to  the  exponent  r(p,  q),  one  obtains  the  two- 
dimensional  D(h,/3 ) singularity  spectrum3 

D(h,0)  = mm(qh+p{l3  + l)-T(p,q))  (7) 

P,Q 

In  practice,  due  to  the  finite  length  and  noise  in  data  sets,  the  accuracy  of  the 
Legendre  transform  might  decrease.  In  order  to  avoid  this  problem,  we  propose  a 
generalization  of  the  canonical  method  of  Chhabra  and  Jensen6  to  two-dimensional 
spectra  which  allows  a direct  computation  of  the  singularity  spectrum  D(h,P). 
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